Advancing chemical safety assessment through an omics-based characterization of the test system-chemical interaction.

Assessing chemical safety is essential to evaluate the potential risks of chemical exposure to human health and the environment. Traditional methods relying on animal testing are being replaced by 3R (reduction, refinement, and replacement) principle-based alternatives, mainly depending on in vitro test methods and the Adverse Outcome Pathway framework. However, these approaches often focus on the properties of the compound, missing the broader chemical-biological interaction perspective. Currently, the lack of comprehensive molecular characterization of the in vitro test system results in limited real-world representation and contextualization of the toxicological effect under study. Leveraging omics data strengthens the understanding of the responses of different biological systems, emphasizing holistic chemical-biological interactions when developing in vitro methods. Here, we discuss the relevance of meticulous test system characterization on two safety assessment relevant scenarios and how omics-based, data-driven approaches can improve the future generation of alternative methods.

Soil Fertility and Phosphorus Leaching in Irrigated Calcareous Soils of the Mediterranean Region.

To ensure soil quality and soil health, it is necessary to improve fertilization practices while minimizing environmental impacts. The aim of this study was to record the state of the art on soil fertility related to fertilization management (organic and/or mineral) and to detect environmental challenges in highly productive fields. A soil survey was set up in a new irrigated area (c. 20 years old), in the north-eastern part of Spain, which is mainly devoted to double annual crop rotations of cereals and maize. The area also supports an important animal rearing activity. The survey covered 733 ha of calcareous soils, owned by 35 farmers. At each farm, fertilization management was recorded, and soil was analyzed for nutrients and heavy metals. Multivariate analyses were performed. Total N, P, Cu and Zn, and available P, Cu, Zn and Mn soil concentrations were associated to the use of organic amendments. Heavy metals concentrations were below established thresholds. Available P (Olsen-P) was identified as an indicator of the previously adopted fertilization management and of the potential of P leaching towards deeper soil layers. Regression analyses were performed. A displacement of available P from the uppermost layer (0-0.3 m) occurs in the breakpoint of 86 mg P kg-1 soil. Preventative actions might be established from 53 mg P kg-1 soil due to the slowdown in P immobilization. Our results reinforce the importance of setting up P threshold soil levels for best practices of fertilization, as a basis for sustainable agriculture intensification.

The in vitro immunomodulatory effect of multi-walled carbon nanotubes by multilayer analysis.

The study of multi-walled carbon nanotube (MWCNT) induced immunotoxicity is crucial for determining hazards posed to human health. MWCNT exposure most commonly occurs via the airways, where macrophages are first line responders. Here we exploit an in vitro assay, measuring dose-dependent secretion of a wide panel of cytokines, as a measure of immunotoxicity following the non-lethal, multi-dose exposure (IC5, IC10 and IC20) to 7 MWCNTs with different intrinsic properties. We find that a tangled structure, and small aspect ratio are key properties predicting MWCNT induced immunotoxicity, mediated predominantly by IL1B cytokine secretion. To assess the mechanism of action giving rise to MWCNT immunotoxicity, transcriptomics analysis was linked to cytokine secretion in a multilayer model established through correlation analysis across exposure concentrations. This reinforced the finding that tangled MWCNTs have greater immunomodulatory potency, displaying enrichment of immune system, signal transduction and pattern recognition associated pathways. Together our results further elucidate how structure, length and aspect ratio, critical intrinsic properties of MWCNTs, are tied to immunotoxicity.

DREAM: an R package for druggability evaluation of human complex diseases.

MOTIVATION: De novo drug development is a long and expensive process that poses significant challenges from the design to the preclinical testing, making the introduction into the market slow and difficult. This limitation paved the way to the development of drug repurposing, which consists in the re-usage of already approved drugs, developed for other therapeutic indications. Although several efforts have been carried out in the last decade in order to achieve clinically relevant drug repurposing predictions, the amount of repurposed drugs that have been employed in actual pharmacological therapies is still limited. On one hand, mechanistic approaches, including profile-based and network-based methods, exploit the wealth of data about drug sensitivity and perturbational profiles as well as disease transcriptomics profiles. On the other hand, chemocentric approaches, including structure-based methods, take into consideration the intrinsic structural properties of the drugs and their molecular targets. The poor integration between mechanistic and chemocentric approaches is one of the main limiting factors behind the poor translatability of drug repurposing predictions into the clinics. RESULTS: In this work, we introduce DREAM, an R package aimed to integrate mechanistic and chemocentric approaches in a unified computational workflow. DREAM is devoted to the druggability evaluation of pathological conditions of interest, leveraging robust drug repurposing predictions. In addition, the user can derive optimized sets of drugs putatively suitable for combination therapy. In order to show the functionalities of the DREAM package, we report a case study on atopic dermatitis. AVAILABILITY AND IMPLEMENTATION: DREAM is freely available at https://github.com/fhaive/dream. The docker image of DREAM is available at: https://hub.docker.com/r/fhaive/dream.

A curated gene and biological system annotation of adverse outcome pathways related to human health.

Adverse outcome pathways (AOPs) are emerging as a central framework in modern toxicology and other fields in biomedicine. They serve as an extension of pathway-based concepts by depicting biological mechanisms as causally linked sequences of key events (KEs) from a molecular initiating event (MIE) to an adverse outcome. AOPs guide the use and development of new approach methodologies (NAMs) aimed at reducing animal experimentation. While AOPs model the systemic mechanisms at various levels of biological organisation, toxicogenomics provides the means to study the molecular mechanisms of chemical exposures. Systematic integration of these two concepts would improve the application of AOP-based knowledge while also supporting the interpretation of complex omics data. Hence, we established this link through rigorous curation of molecular annotations for the KEs of human relevant AOPs. We further expanded and consolidated the annotations of the biological context of KEs. These curated annotations pave the way to embed AOPs in molecular data interpretation, facilitating the emergence of new knowledge in biomedicine.

ESPERANTO: a GLP-field sEmi-SuPERvised toxicogenomics metadAta curatioN TOol.

SUMMARY: Biological data repositories are an invaluable source of publicly available research evidence. Unfortunately, the lack of convergence of the scientific community on a common metadata annotation strategy has resulted in large amounts of data with low FAIRness (Findable, Accessible, Interoperable and Reusable). The possibility of generating high-quality insights from their integration relies on data curation, which is typically an error-prone process while also being expensive in terms of time and human labour. Here, we present ESPERANTO, an innovative framework that enables a standardized semi-supervised harmonization and integration of toxicogenomics metadata and increases their FAIRness in a Good Laboratory Practice-compliant fashion. The harmonization across metadata is guaranteed with the definition of an ad hoc vocabulary. The tool interface is designed to support the user in metadata harmonization in a user-friendly manner, regardless of the background and the type of expertise. AVAILABILITY AND IMPLEMENTATION: ESPERANTO and its user manual are freely available for academic purposes at https://github.com/fhaive/esperanto. The input and the results showcased in Supplementary File S1 are available at the same link.

Oribatid mites in different Mediterranean crop rotations fertilized with animal droppings.

Unsustainable soil management is one of the drivers of soil degradation, but impact assessment requires the development of indicators. Oribatids might be considered as early indicators of disturbances due to the stability of their community. The aim of this study was to investigate the feasibility of oribatids as bioindicators of sustainable agricultural practices. Under a dry Mediterranean climate, three fertilization experiments - two under a two-crop rotation system and one under maize monoculture and established 12 years earlier - were sampled 3× for oribatid identification during the last annual cropping cycle. The hypothesis was that different nutrient and crop managements affect the number of oribatid species and individuals present, and these parameters could be used as indicators of soil degradation. In total, 18 oribatid species were identified, and 1974 adult individuals were recovered. Maximum abundance was found prior to sowing. Pig slurry (PS) vs. control, and dairy cattle manure (CM) vs. mineral fertilization increased oribatid abundance. This increase was evident when the average applied rates with PS were ca. 2 Mg of organic matter (OM) ha- 1 yr- 1, or higher than ca. 4 Mg OM ha- 1 yr- 1 for CM. When the preceding crop was wheat and PS or CM were used, Oribatula (Zygoribatula) excavata (which reproduces sexually) predominated. In maize monoculture fertilized with CM, Tectocepheus sarekensis and Acrotritia ardua americana (which can reproduce through parthenogenesis) prevailed vs. Oribatula, which indicated a heavily disturbed soil. Under this specific Mediterranean environment, the predominance of certain parthenogenic oribatid species and the number of individuals provide advanced warning on soil degradation.

KNeMAP: a network mapping approach for knowledge-driven comparison of transcriptomic profiles.

MOTIVATION: Transcriptomic data can be used to describe the mechanism of action (MOA) of a chemical compound. However, omics data tend to be complex and prone to noise, making the comparison of different datasets challenging. Often, transcriptomic profiles are compared at the level of individual gene expression values, or sets of differentially expressed genes. Such approaches can suffer from underlying technical and biological variance, such as the biological system exposed on or the machine/method used to measure gene expression data, technical errors and further neglect the relationships between the genes. We propose a network mapping approach for knowledge-driven comparison of transcriptomic profiles (KNeMAP), which combines genes into similarity groups based on multiple levels of prior information, hence adding a higher-level view onto the individual gene view. When comparing KNeMAP with fold change (expression) based and deregulated gene set-based methods, KNeMAP was able to group compounds with higher accuracy with respect to prior information as well as is less prone to noise corrupted data. RESULT: We applied KNeMAP to analyze the Connectivity Map dataset, where the gene expression changes of three cell lines were analyzed after treatment with 676 drugs as well as the Fortino et al. dataset where two cell lines with 31 nanomaterials were analyzed. Although the expression profiles across the biological systems are highly different, KNeMAP was able to identify sets of compounds that induce similar molecular responses when exposed on the same biological system. AVAILABILITY AND IMPLEMENTATION: Relevant data and the KNeMAP function is available at: https://github.com/fhaive/KNeMAP and 10.5281/zenodo.7334711.

Toxicogenomics Data for Chemical Safety Assessment and Development of New Approach Methodologies: An Adverse Outcome Pathway-Based Approach.

Mechanistic toxicology provides a powerful approach to inform on the safety of chemicals and the development of safe-by-design compounds. Although toxicogenomics supports mechanistic evaluation of chemical exposures, its implementation into the regulatory framework is hindered by uncertainties in the analysis and interpretation of such data. The use of mechanistic evidence through the adverse outcome pathway (AOP) concept is promoted for the development of new approach methodologies (NAMs) that can reduce animal experimentation. However, to unleash the full potential of AOPs and build confidence into toxicogenomics, robust associations between AOPs and patterns of molecular alteration need to be established. Systematic curation of molecular events to AOPs will create the much-needed link between toxicogenomics and systemic mechanisms depicted by the AOPs. This, in turn, will introduce novel ways of benefitting from the AOPs, including predictive models and targeted assays, while also reducing the need for multiple testing strategies. Hence, a multi-step strategy to annotate AOPs is developed, and the resulting associations are applied to successfully highlight relevant adverse outcomes for chemical exposures with strong in vitro and in vivo convergence, supporting chemical grouping and other data-driven approaches. Finally, a panel of AOP-derived in vitro biomarkers for pulmonary fibrosis (PF) is identified and experimentally validated.

The potential of a data centred approach & knowledge graph data representation in chemical safety and drug design.

Big Data pervades nearly all areas of life sciences, yet the analysis of large integrated data sets remains a major challenge. Moreover, the field of life sciences is highly fragmented and, consequently, so is its data, knowledge, and standards. This, in turn, makes integrated data analysis and knowledge gathering across sub-fields a demanding task. At the same time, the integration of various research angles and data types is crucial for modelling the complexity of organisms and biological processes in a holistic manner. This is especially valid in the context of drug development and chemical safety assessment where computational methods can provide solutions for the urgent need of fast, effective, and sustainable approaches. At the same time, such computational methods require the development of methodologies suitable for an integrated and data centred Big Data view. Here we discuss Knowledge Graphs (KG) as a solution to a data centred analysis approach for drug and chemical development and safety assessment. KGs are knowledge bases, data analysis engines, and knowledge discovery systems all in one, allowing them to be used from simple data retrieval, over meta-analysis to complex predictive and knowledge discovery systems. Therefore, KGs have immense potential to advance the data centred approach, the re-usability, and informativity of data. Furthermore, they can improve the power of analysis, and the complexity of modelled processes, all while providing knowledge in a natively human understandable network data model.

Representing and describing nanomaterials in predictive nanoinformatics.

Engineered nanomaterials (ENMs) enable new and enhanced products and devices in which matter can be controlled at a near-atomic scale (in the range of 1 to 100 nm). However, the unique nanoscale properties that make ENMs attractive may result in as yet poorly known risks to human health and the environment. Thus, new ENMs should be designed in line with the idea of safe-and-sustainable-by-design (SSbD). The biological activity of ENMs is closely related to their physicochemical characteristics, changes in these characteristics may therefore cause changes in the ENMs activity. In this sense, a set of physicochemical characteristics (for example, chemical composition, crystal structure, size, shape, surface structure) creates a unique 'representation' of a given ENM. The usability of these characteristics or nanomaterial descriptors (nanodescriptors) in nanoinformatics methods such as quantitative structure-activity/property relationship (QSAR/QSPR) models, provides exciting opportunities to optimize ENMs at the design stage by improving their functionality and minimizing unforeseen health/environmental hazards. A computational screening of possible versions of novel ENMs would return optimal nanostructures and manage ('design out') hazardous features at the earliest possible manufacturing step. Safe adoption of ENMs on a vast scale will depend on the successful integration of the entire bulk of nanodescriptors extracted experimentally with data from theoretical and computational models. This Review discusses directions for developing appropriate nanomaterial representations and related nanodescriptors to enhance the reliability of computational modelling utilized in designing safer and more sustainable ENMs.

Biomarkers of nanomaterials hazard from multi-layer data.

There is an urgent need to apply effective, data-driven approaches to reliably predict engineered nanomaterial (ENM) toxicity. Here we introduce a predictive computational framework based on the molecular and phenotypic effects of a large panel of ENMs across multiple in vitro and in vivo models. Our methodology allows for the grouping of ENMs based on multi-omics approaches combined with robust toxicity tests. Importantly, we identify mRNA-based toxicity markers and extensively replicate them in multiple independent datasets. We find that models based on combinations of omics-derived features and material intrinsic properties display significantly improved predictive accuracy as compared to physicochemical properties alone.

Characterization of ENM Dynamic Dose-Dependent MOA in Lung with Respect to Immune Cells Infiltration.

The molecular effects of exposures to engineered nanomaterials (ENMs) are still largely unknown. In classical inhalation toxicology, cell composition of bronchoalveolar lavage (BAL) is a toxicity indicator at the lung tissue level that can aid in interpreting pulmonary histological changes. Toxicogenomic approaches help characterize the mechanism of action (MOA) of ENMs by investigating the differentially expressed genes (DEG). However, dissecting which molecular mechanisms and events are directly induced by the exposure is not straightforward. It is now generally accepted that direct effects follow a monotonic dose-dependent pattern. Here, we applied an integrated modeling approach to study the MOA of four ENMs by retrieving the DEGs that also show a dynamic dose-dependent profile (dddtMOA). We further combined the information of the dddtMOA with the dose dependency of four immune cell populations derived from BAL counts. The dddtMOA analysis highlighted the specific adaptation pattern to each ENM. Furthermore, it revealed the distinct effect of the ENM physicochemical properties on the induced immune response. Finally, we report three genes dose-dependent in all the exposures and correlated with immune deregulation in the lung. The characterization of dddtMOA for ENM exposures, both for apical endpoints and molecular responses, can further promote toxicogenomic approaches in a regulatory context.

Influence of FLG loss-of-function mutations in host-microbe interactions during atopic skin inflammation.

BACKGROUND: Loss-of-function mutations in the filaggrin (FLG) gene directly alter skin barrier function and critically influence atopic inflammation. While skin barrier dysfunction, Th2-associated inflammation and bacterial dysbiosis are well-known characteristics of atopic dermatitis (AD), the mechanisms interconnecting genotype, transcriptome and microbiome remain largely elusive. OBJECTIVE: In-depth analysis of FLG genotype-associated skin gene expression alterations and host-microbe interactions in AD. METHODS: Multi-omics characterization of a cohort of AD patients carrying heterozygous loss-of-function mutations in the FLG gene (ADMut) (n = 15), along with matched wild-type (ADWt) patients and healthy controls. Detailed clinical characterization, microarray gene expression and 16 S rRNA-based microbial marker gene data were generated and analyzed. RESULTS: In the context of filaggrin dysfunction, the transcriptome was characterized by dysregulation of barrier function and water homeostasis, while the lesional skin of ADWt demonstrated the specific upregulation of pro-inflammatory cytokines and T-cell proliferation. S. aureus dominated the microbiome in both patient groups, however, shifting microbial communities could be observed when comparing healthy with non-lesional ADWt or ADMut skin, offering the opportunity to identify microbe-associated transcriptomic signatures. Moreover, an AD core signature with 28 genes, including CCL13, CCL18, BTC, SCIN, RAB31 and PCLO was identified. CONCLUSIONS: Our integrative approach provides molecular insights for the concept that FLG loss-of-function mutations are a genetic shortcut to atopic inflammation and unravels the complex interplay between genotype, transcriptome and microbiome in the human holobiont.

Seven years of pig slurry fertilization: impacts on soil chemical properties and the element content of winter barley plants.

Intensive pig farming produces large amounts of slurry, which is applied to agricultural soils as fertilizer. A 7-year field study was performed to check the effect of pig slurry on soil properties and on the accumulation of some essential nutrients and heavy metals in a calcareous silty-loam soil (0-0.3 m) and in barley (Hordeum vulgare L.) plants in two cropping seasons with contrasting amounts of rainfall. Five fertilization treatments, control (no N applied), mineral fertilizer (90 kg N ha-1), and different N doses of pig slurry (146, 281, 534 kg N ha-1), were applied at sowing of a barley crop. Organic carbon, available P and K, and total P in soil increased with slurry dose. No differences were found in Co, Cr, Fe, Mn, Ni, and Pb soil concentrations. Slurries increased Cu, Mn, and Zn extractions and plant concentrations of P in straw and Zn in grain. However, the lowest slurry rate was able to maintain the highest grain yields while improving fertility. The results of this research study support the sustainability of pig slurry fertilization at appropriate rates in relation to soil chemical quality.

Integrated Network Pharmacology Approach for Drug Combination Discovery: A Multi-Cancer Case Study.

Despite remarkable efforts of computational and predictive pharmacology to improve therapeutic strategies for complex diseases, only in a few cases have the predictions been eventually employed in the clinics. One of the reasons behind this drawback is that current predictive approaches are based only on the integration of molecular perturbation of a certain disease with drug sensitivity signatures, neglecting intrinsic properties of the drugs. Here we integrate mechanistic and chemocentric approaches to drug repositioning by developing an innovative network pharmacology strategy. We developed a multilayer network-based computational framework integrating perturbational signatures of the disease as well as intrinsic characteristics of the drugs, such as their mechanism of action and chemical structure. We present five case studies carried out on public data from The Cancer Genome Atlas, including invasive breast cancer, colon adenocarcinoma, lung squamous cell carcinoma, hepatocellular carcinoma and prostate adenocarcinoma. Our results highlight paclitaxel as a suitable drug for combination therapy for many of the considered cancer types. In addition, several non-cancer-related genes representing unusual drug targets were identified as potential candidates for pharmacological treatment of cancer.

Nextcast: A software suite to analyse and model toxicogenomics data.

The recent advancements in toxicogenomics have led to the availability of large omics data sets, representing the starting point for studying the exposure mechanism of action and identifying candidate biomarkers for toxicity prediction. The current lack of standard methods in data generation and analysis hampers the full exploitation of toxicogenomics-based evidence in regulatory risk assessment. Moreover, the pipelines for the preprocessing and downstream analyses of toxicogenomic data sets can be quite challenging to implement. During the years, we have developed a number of software packages to address specific questions related to multiple steps of toxicogenomics data analysis and modelling. In this review we present the Nextcast software collection and discuss how its individual tools can be combined into efficient pipelines to answer specific biological questions. Nextcast components are of great support to the scientific community for analysing and interpreting large data sets for the toxicity evaluation of compounds in an unbiased, straightforward, and reliable manner. The Nextcast software suite is available at: ( https://github.com/fhaive/nextcast).

Microarray Data Preprocessing: From Experimental Design to Differential Analysis.

DNA microarray data preprocessing is of utmost importance in the analytical path starting from the experimental design and leading to a reliable biological interpretation. In fact, when all relevant aspects regarding the experimental plan have been considered, the following steps from data quality check to differential analysis will lead to robust, trustworthy results. In this chapter, all the relevant aspects and considerations about microarray preprocessing will be discussed. Preprocessing steps are organized in an orderly manner, from experimental design to quality check and batch effect removal, including the most common visualization methods. Furthermore, we will discuss data representation and differential testing methods with a focus on the most common microarray technologies, such as gene expression and DNA methylation.

Supervised Methods for Biomarker Detection from Microarray Experiments.

Biomarkers are valuable indicators of the state of a biological system. Microarray technology has been extensively used to identify biomarkers and build computational predictive models for disease prognosis, drug sensitivity and toxicity evaluations. Activation biomarkers can be used to understand the underlying signaling cascades, mechanisms of action and biological cross talk. Biomarker detection from microarray data requires several considerations both from the biological and computational points of view. In this chapter, we describe the main methodology used in biomarkers discovery and predictive modeling and we address some of the related challenges. Moreover, we discuss biomarker validation and give some insights into multiomics strategies for biomarker detection.

Unsupervised Algorithms for Microarray Sample Stratification.

The amount of data made available by microarrays gives researchers the opportunity to delve into the complexity of biological systems. However, the noisy and extremely high-dimensional nature of this kind of data poses significant challenges. Microarrays allow for the parallel measurement of thousands of molecular objects spanning different layers of interactions. In order to be able to discover hidden patterns, the most disparate analytical techniques have been proposed. Here, we describe the basic methodologies to approach the analysis of microarray datasets that focus on the task of (sub)group discovery.